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THE QUILLENSUSLIN PACKAGE FOR MACAULAY2 



BRETT BARWICK AND BRANDEN STONE 



Abstract. The QuillenSuslin package for Macaulay2 provides the abiUty 
to compute a free basis for a projective module over a polynomial ring with 
coefficients in Q, Z, or Z/p for a prime integer p. A brief description of the 
underlying algorithm and the related tools are given. 



1. Introduction 

\^ ' In 1955, J. P. Serre posed the following question: Do there exist finitely generated 

projective fc[a;i, . . . ,a;„] modules, with k a field, which are not free? [Ser55] This 
question was known as "Serre's Problem" and the question in its full generality 
remained open for 21 years until it was resolved independently by D. Quillen and 
2 ' A. A. Suslin in 1976, resulting in the following well-known theorem. 

Theorem 1 (Quillen-Suslin, 1976 [Qui76, Sus76]). Let S = k[xi, . . . ,a;„], with k a 
field. Then every finitely generated projective S-module is free. 

^ ' However, the proofs given were not entirely constructive, and it was not until 

Cn ! the early 1990's that papers such as [FG90, LS92, LWOO] began giving fully con- 

00 ' structive versions of the proof. In 1992, A. Logar and B. Sturmfels [LS92] published 

the algorithmic proof of the Quillen-Suslin Theorem that forms the basis for the 

methods in QuillenSuslin. In their paper, Logar and Sturmfels describe, via the 

^-^ I technique of completion of unimodular rows, how to construct a free generating 

set for a projective module over C[a;i, . . . ,x„]. One can extend these constructive 
techniques to work over more general coefficient rings such as Q, Z, and Z/pZ, for p 
a prime integer. Descriptions of some of these more general techniques can be found 
in [Fab09, Ch. 2], and we have implemented these algorithms, with some modifi- 
j^ ' cations, in our QuillenSuslin package for Macaulay2 [GS]. In the next section we 

JH I will give some preliminary definitions and results which reduce the statement of 

the Quillen-Suslin Theorem to a more concrete matrix theoretic problem concern- 
ing the completion of unimodular rows over polynomial rings to square invertible 
matrices. 

2. Preliminaries 

In this section, R will denote a commutative ring and AI will denote a finitely 
generated i?-module. We say that M is a projective i?-module if it is a direct 
summand of a free module. Similarly, we define a slightly stronger notion by saying 
that M is stably free if there exists some to > such that M © _R'" is free. A 
module M is stably free if and only if it is isomorphic to the kernel of a surjective 
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i?-linear map </> : i?" — > i?™ for some m < n. Since <j> surjects onto a free module, 
we know that this map spHts and admits a right inverse ip : i?™ — > i?" so that 
(pip = idijTTi . Therefore a matrix representing (p is right invertible, and we caU 
such a right invertible matrix over R unimodular. Using this terminology, it is not 
difficult to show that the Quillen-Suslin Theorem as stated above is equivalent to 
the following matrix theoretic statement about unimodular matrices. 

Theorem 2 (Quillen-Suslin, restatement [LS92, Theorem 1.1]). Let S = R[xi, . . . , x„ 
with R a principal ideal domain and let U £ Matmxn{R) be a unimodular matrix 
over S with m < n. Then there exists a unimodular matrix V £ MatnxniS) such 
that 



UV = 



A matrix V which satisfies the properties in the above theorem is said to solve the 
unimodular matrix problem for U. Given such a V, the first m rows of the invertible 
matrix V~^ are the same as the original matrix U. Therefore we see that proving 
the Quillen-Suslin Theorem is equivalent to showing that any unimodular matrix 
can be completed to a square invertible matrix over the polynomial ring. As it 
turns out, it suffices to show that the unimodular row problem can be solved, that is. 
Theorem 2 holds for unimodular row vectors [LS92]. For more details concerning 
this equivalent formulation of the Quillen-Suslin Theorem, we refer the interested 
reader to the excellent book of Lam [Lam06] . 

3. The Logar-Sturmfels Algorithm 

Before describing the general algorithm, we mention that QuillenSuslin contains 
several shortcut methods, as described in [Fab09, Sect. 2.2] which allow us to 
quickly solve the unimodular row problem for a row satisfying certain properties, 
often allowing us to avoid the worst-case general algorithm. These shortcut methods 
are automatically used as soon as they are applicable at any point during the 
methods computeFreeBasis, completeMatrix, and qsAlgorithm. 

The main idea behind the Logar-Sturmfels algorithm is to iteratively reduce the 
number of variables involved by one in a unimodular row /, obtaining a unimodular 
row / over the coefficient ring, which is a PID. We can then use a simple algorithm 
based on the Smith normal form of / to construct a final unimodular matrix U so 
that fU =[10 • • • 0] . Multiplying together all of the matrices used during the 
process, one can construct a unimodular matrix which solves the unimodular row 
problem for the original row /. 

The process of eliminating a variable from a unimodular matrix is organized 
into three main steps: the normalization step, the "local loop", and the patching 
step. Below is a brief description of each step, as well as a demonstration of the 
corresponding commands in the QuillenSuslin package. We will work over the poly- 
nomial ring S ~ 1i[x,y\ and consider the unimodular row f = [x^ , 2y-\-\,x^y^ + y\ 
over iS*. 

il : loadPackage "QuillenSuslin"; 

12 : S = ZZ[x,y] ; 
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13 : f = matrix {{x~2,2*y+l ,x~5*y-2+y» 
o3 = I x2 2y+l x5y2+y I 
1 3 

03 : Matrix S < S 

We can use the command isUnimodular to check that this row is indeed uni- 
modular over S. 

i4 : isUnimodular f 

04 = true 

In order to ehminate the variable y from this unimodular row, we will construct 
a unimodular matrix U so that multiplying f ^U is the same as evaluating f when 
2/ = 0. We first demonstrate the normalization step. 

3.1. Normalization Step. Since Horrocks' Theorem (see Theorem 4 below) re- 
quires a monic polynomial, we must first construct a unimodular matrix U and an 
invertible change of variables xi <-> Xi so that the first entry of f *U is monic in Xn- 
The normalization step is based on the following result, which has been slightly 
restated for our purposes. 

Lemma 3 ([VS76, Lemma 10.6]). Let R be a Noetherian ring, n and m natural 
numbers, S = R[xi, . . . , Xn], m > dimi? + 2, and f = [/i /2 • • ■ fm] o, unimod- 
ular row over S . Then there exists an m x m unimodular matrix U over S, and 
an invertible change of variables Xi <-> Xi, so that after applying the change of 
variables the first entry of fU is monic in Xn when viewed as a polynomial in 
R[Xi, . . . , X„^i][X„]. 

A constructive version of this result is implemented in the method changeVar, 
and is used as in the following example. 

i5 : (Ul ,subs , invSubs) = changeVar (f ,{x,y}-) 
o5=(|100|,|yx|,|yx|) 

I 1 I 

I 1 I 

05 : Sequence 

i6 : f = sub (f*Ul, subs) 

06 = I y2 2x+l x2y5+x I 

1 3 

o6 : Matrix S < S 

Notice that since the first entry of the row f was already monic in a;, the method 
simply returned a permutation of the variables interchanging x and y so that the 
first entry of the new row would be monic in y. Now that the first entry of the row 
is monic in the variable we are trying to eliminate, we may proceed to the "local 
loop." 

3.2. Local Loop. The purpose of the local loop is to compute a collection of local 
solutions to the unimodular row problem for f . The local loop is based on the 
following result of Horrocks. 
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Theorem 4 (Horrocks, [Rot09, Prop. 4.98]). Consider the polynomial ring B[y], 
where B is a local ring, and let f ~ [fi f2 ■ • ■ fm] be a unimodular row over B[y]. 
If some fi is monic in y, then there exists a unimodular mx m matrix U over B[y] 
so that fU =[10 ■■■ 0] . 

In order to eliminate the last variable Xn in a unimodular row over a polynomial 
ring R[xi, . . . ,x„], the local loop proceeds in the following way: First set / = (0) 
in A = R[xi, . . . , Xn-i]- Now while I ^ A, the i*'' iteration of the loop is 

(1) Find a maximal ideal trii in A containing /. 

(2) Apply Horrocks' Theorem to the row /, viewed as a unimodular row over 
Ami i^n], to find a unimodular matrix Li over Ajm [xn] which solves the uni- 
modular row problem for / (we call this L,; a local solution to the unimodular 
row problem for /). 

(3) Let di denote the common denominator for all of the elements in the matrix 

(4) Set I = I+{di). 

li I y^ A, then we repeat the loop. Otherwise we are able to stop and go on to the 
patching step. Notice that since we are creating a strictly ascending chain of ideals 
(di) C {di,d2) C • • • in the Noetherian ring A, this loop must terminate in a finite 
number of steps with (di, . . . ,dk) = A for some integer k. 

In our example, we use the method getMaxIdeal to first find an arbitrary max- 
imal ideal in Z[a;], and we set mi = (2,a;). Using the method horrocks, we can 
compute a unimodular matrix LI over (Z[a;](2,a;)) [y] so that f *L1 = I 10 I . 

17 : A = ZZ[x] ; 

18 : ml = getMaxIdeal (Ideal (0_A) ,{x}) 
o8 = Ideal (2, x) 

08 : Ideal of A 

19 : LI = horrocks (f ,y,ml) 

09 = I 1 I 

I l/(2x+l) -y2/(2x+l) (-x2y5-x)/(2x+l) I 
10 1 I 

3 3 

09 : Matrix (frac S) < (frac S) 

Since di = 2a;+l is a common denominator for the entries of Ll and (2.T+1) ^ Z[x], 
we use getMaxIdeal again to find a maximal ideal containing 2a; + 1, and we set 
m2 = (3,x — 1). We use horrocks a second time to compute a new local solution 
L2 with common denominator d2 = x. 

110 : m2 = getMaxIdeal(sub(ldeal(2*x+l) ,A) ,{x}) 

010 = Ideal (x - 1, 3) 
olO : Ideal of A 

111 : L2 = horrocks (f,y,m2) 
oil = I -xy3 xy5+l 2x2y3+xy3 I 

10 1 I 

I 1/x -y2/x (-2x-l)/x I 
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3 3 

oil : Matrix (frac S) < (frac S) 

il2 : sub(ideal(2*x+l,x) ,S) == ideal(l_S) 

012 = true 

Since {di,d2) = (2x + l,a;) = Z[a;], we are able to exit the local loop and proceed 
to the patching step. 

3.3. Patching Step. Loosely speaking, the patching step involves multiplying 
slight variations of the local solutions Li, . . . ,Lk together in a clever way so that the 
product [/ is a unimodular matrix over the polynomial ring R[xi, . . . , a;„] and mul- 
tiplying / times U is equivalent to evaluating / when Xn = 0, thereby eliminating 
one of the variables in the row /. 

Following along with our example, we use the method patch applied to our list 
{L1,L2} of local solutions and we specify that y is the variable that we want to 
eliminate. 

il3 : U = patch({Ll,L2},y) 

013 = I -32x6y5+l 8x5y3 I 

I 16x5y7-8x4y7+4x3y7+2xy2-y2 1 -4x4y5+2x3y5-x2y5 I 
I -4xy2 1 I 

3 3 

013 : Matrix S < S 

il4 : f*U 

oi4 = I 2x+l X I 

1 3 

014 : Matrix S < S 

We can see that multiplying the row f times the unimodular matrix U is equivalent 
to evaluating f when y = (keeping in mind that the variables x and y were 
interchanged during the normalization step). 

4. Core methods in the QuillenSuslin package 

The method qsAlgorithm automates all of the above computations for com- 
puting a solution to the unimodular matrix problem, and automatically applies the 
shortcut methods in [Fab09] when possible. We demonstrate the use of qsAlgorithm 
by finding a solution to the unimodular row problem for the row f = [y'^ , 2x + 
1 , x^y^ + x] given earlier. 

il5 : U = qsAlgorithm f 

015 = I 2x3y2 -2x5y2+l -2x8y4-2x3y3 I 

I 1 -x2 -x5y2-y I 

I -2 2x2 2x5y2+2y+l I 

3 3 

015 : Matrix S < S 

il6 : det U 

016 = -1 
ol6 = S 
il7 : f*U 
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0l7 = I 1 I 

1 3 

017 : Matrix S < S 

The package also contains a method completeMatrix which completes a uni- 
modular matrix over a polynomial ring to a square invertible matrix. Again we 
demonstrate its use on the unimodular row f . 
il8 : C = completeMatrix f 

018 = {0} I y2 2x+l x2y5+x I 

{2} I 1 -2x2y3 I 

{7} I 2 1 I 

3 3 

018 : Matrix S < S 

il9 : det C 

019 = -1 

019 : S 

Finally, we give an example to demonstrate the method computeFreeBasis, 
which computes a free generating set for a projective module. We define K = 
kerf, which we can check is a projective Z[a;, j/]-module by using the command 
isProjective. 

i20 : K = ker f 

020 = image I 2x+l 2x3y3+x2y3 2x2y5-l -x2y5-x I 

I -y2 X y2 I 
I -2x-l -2y2 y2 I 
3 

020 : S-module, submodule of S 

i21 : isProjective K 

021 = true 

i22 : mingens K 

022 = I -2x-l 2x3y3+x2y3 x2y5+x I 

I y2 X I 
I -2x-l -y2 I 
3 3 

022 : Matrix S < S 

i23 : B = computeFreeBasis K 

023 = I -2x2y5+l 2x3y3+x2y3 I 

I -y2 X I 

I 2y2 -2x-l I 

3 2 

023 : Matrix S < S 

i24 : syz B 

024 = 

2 
o24 : Matrix S < 
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125 : image B == K 

o25 = true 
From the Macaulay2 output, we can see that the native command mingens does 
not produce a free generating set for K. but computeFreeBasis does produce a set 
of 2 generators for K with no relations, demonstrating that K is free. 
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